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ABSTRACT 

We present the results of a survey for super Lyman limit systems (SLLS; denned to be absorbers 
with 19.0 <log Nhi < 20.3 cm -2 ) from a large sample of high resolution spectra acquired using the 
Keck and Magellan telescopes. Specifically, we present 47 new SLLS from 113 QSO sightlines. We focus 
on the neutral hydrogen frequency distribution /hi (A, X) of the SLLS and its moments, and compare 
these results with the Lyman-a forest and the damped Lyman alpha systems (DLA; absorbers with 
logiVf27> 20.3 cm -2 ). We find that that fm(N,X) of the SLLS can be reasonably described with a 



power-law of index a = — 1.43_q i6 or a = ~^--^-o.2i depending on whether we set the lower Ahi bound 
for the analysis at 10 cm -2 or 10 193 cm -2 , respectively. The results indicate a flattening in the slope 
of fm(N,X) between the SLLS and DLA. We find little evidence for redshift evolution in the shape 
of fm(N,X) for the SLLS over the redshift range of the sample 1.68 < z < 4.47 and only tentative 
evidence for evolution in the zeroth moment of /hi (N, X), the line density ^slls(A). We introduce the 
observable distribution function 0(N, X) and its moment, which elucidates comparisons of H I absorbers 
from the Lyman-a forest through to the DLA. We find that a simple three parameter function can fit 
0(N, X) over the range 17.0 <\ogN HI < 22.0. We use these results to predict that f m (N,X) must 
show two additional inflections below the SLLS regime to match the observed fm(N, X) distribution 
of the Lyman-a forest. Finally, we demonstrate that SLLS contribute a minor fraction (s=s 15%) of the 
universe's hydrogen atoms and, therefore, an even smaller fraction of the mass in predominantly neutral 
gas. 

Subject headings: quasars: absorption lines - intergalactic medium 

1. INTRODUCTION 



For nearly three decades, the study of absorption line 
systems towards distant quasars has addressed a wide range 
of astrophysical and cosmological issues. These systems 
are typically classified according to their neutral hydro- 
gen content: the Lya forest absorbers with log Nhi< 17.2 
cm -2 , the Lyman limit systems (LLS) with 17.2 <log Nhi< 
20.3 cm~ 2 , and the damped Lyman alpha systems (DLA) 
with \ogN HI > 20.3 cm -2 . Both the Lya forest and DLA 
absorbers have received considerable attention. The Lya 
forest can be used to constrain cosmological parameters 
through a number of methods, such as through studies 
of the flux power spectrum (Croft et al. 2002; McDon- 
ald et al. 2005), the mean flux decrement (Tytler et al. 
2004), or the distribution in column density and veloc- 
ity width of the absorbers (e.g. Kirkman & Tytler 1997; 
Kim et al. 2002). The DLA trace the bulk of the neutral 
gas at high redshift and are believed to be the progeni- 
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tors of modern day galaxies (Wolfe, Gawiser & Prochaska 
2005). Large statistical samples of both the Lya forest 
and the DLA are readily observed because their Ahi is 
easily determined. For the Lya forest, the JVhi is deter- 
mined directly from Voigt profile fits to the Lya line which 
is dominated by the Maxwellian profile, or from higher or- 
der Lyman transitions when Lya is saturated. For the 
DLA, the Lorentzian component of the Voigt profile gives 
pronounced damping wings, allowing for accurate TVhi de- 
terminations from spectra even at low signal to noise or 
resolution. 

Recently, Prochaska, Herbert-Fort, & Wolfe (2005, here- 
after PHW2005) analyzed the thousands of spectra from 
the Sloan Digital Sky Survey (SDSS) Data Release 3, and 
determined the Ahi frequency distribution for over 500 
DLA systems. In contrast, there exists comparatively lit- 
tle study of the LLS. Surveys for LLS absorption (Tytler 
1982; Sargent, Steidel & Boksenberg 1989; Lanzetta 1991; 
Storrie-Lombardi et al. 1994; Stengler-Larrea et al. 1995) 
have concentrated primarily on the frequency of absorp- 
tion with redshift (frequently expressed as dn/dz or dN/dz, 
but we adopt the notation l(z) for the line density which 
was introduced by PHW2005), but not on the Nm value 
of the LLS. These surveys generally included the full range 
of H I column density of log Nhi > 17.2 cm -2 in £(z). As 
such, the surveys also contained the £ (z) of DLA systems. 

Our ignorance of LLS largely stems from the fact that 
accurate determinations of Ahi for the LLS are difficult 
compared to the Lya forest and the DLA. In part, this is 
because the LLS systems represent the flat portion of the 
curve of growth for the Lya transition and a precise Ahi 
measurement requires high resolution observations (e.g. 
Steidel 1990). For most LLS, the Ahi is determined by 
using the information from the Lyman limit, cither by 
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looking at the differential flux level above and below the 
limit, or by the spacing of the lines as they approach the 
Lyman limit, or both (e.g. Buries & Tytler 1998). The dif- 
ferential flux method requires a detailed understanding of 
the continuum flux level which is difficult to obtain in high 
rcdshift QSOs, and the line spacing technique requires a 
precise model for the hydrogen velocity structure, which 
can generally only be inferred from the metal lines asso- 
ciated with the LLS (e.g. Kirkman et al. 2003). These 
challenges are particularly acute for low-resolution obser- 
vations where the Lyman series is poorly resolved. 

Because the frequency of intersecting a LLS is observed 
to be of order one per sightline at redshifts z > 2.5 (Storrie- 
Lombardi et al. f994), a relatively large QSO sample is re- 
quired in comparison with Lya forest studies. At redshifts 
z > 3.25, the SDSS provides just such a sample (Adelman- 
McCarthy et al. 2006), offering many thousands of sight - 
lincs suitable for analysis of LLS. Unfortunately, the SDSS 
spectra are of too poor a resolution to provide useful con- 
straints on the Ahi for the LLS which are optically thick, 
i.e. log Ay/ > 17.5 cm~ 2 . To address this nearly three 
decade gap in Ahi, we have pursued a high-resolution sur- 
vey for LLS using telescopes in both the northern and 
southern hemisphere. The goal of this survey is to obtain 
high-resolution, high-SNR spectra of at least 100 LLS, 
with full coverage over the LLS Ahi range. In this pa- 
per, we discuss a sub-sample of this survey, namely the 
LLS with \ogN HI > 19.0 cm" 2 . The lower bound on A ffl 
is chosen such that systems with this column density can 
be easily identified and analyzed in moderate to high res- 
olution data (FWHM < 50 km s _1 ; Dessauges-Zavadsky 
et al. 2003). The Ahi values for these LLS are deter- 
mined from the damping wings present in the Lya line, in 
a fashion analogous with DLA analysis in lower resolution 
spectra. Previous work on this subset of the LLS popu- 
lation referred to the absorption systems as 'sub-DLAs'. 
Because the majority of these absorbers are likely to be 
predominantly ionized (e.g. Viegas 1995; Prochaska 1999), 
we adopt the nomenclature of Prochaska & Herbert-Fort 
(2004) for those LLS exhibiting \ogN HI > 19.0 cm" 2 : the 
'super Lyman limit' or SLLS absorbers. 

The fundamental measure of a class of QAL systems is 
the H I column density frequency distribution /hi (A, A), 
defined to be the number of absorbers in the column den- 
sity interval (A, A + dN) identified along the cosmological 
distance interval (A, A + dX) with dX = t#jy(1 + zfdz. 
This quantity, the absorption distance, is generally evalu- 
ated across a redshift interval Az. The /hi(A, A) distribu- 
tion for QAL systems is analogous to the luminosity func- 
tion of galaxy surveys. Moments of /hi (A, A) give impor- 
tant quantities such as the line density of absorption sys- 
tems, and the mass density of H I atoms. The /hi (A, A) 
distribution is the starting point for assessing the baryonic 
mass density of these absorbers as well as their cosmolog- 
ical metal budget. The frequency distribution of LLS is 
of particular interest because the Ahi = 10 17 — 10 20 cm~ 2 
interval is expected to include the Ahi value where QAL 
systems transition from primarily neutral gas to predom- 
inantly ionized gas (e.g. Zheng & Miralda-Escude 2002, ; 
PHW05). Furthermore, Prochaska et al. (2006) have ar- 
gued that the LLS can constitute a considerable fraction 
of the metals in the young universe (see also Peroux et al. 



2006). At z ~ 2, a census for metals which includes the 
DLA, stars in high- z galaxies, and the IGM falls short by 
up to 70% of the total predicted metal mass density (e.g. 
Bouche, Lehnert & Peroux 2006). 

The primary goal for this study of the SLLS is to ex- 
tend the statistics we use to describe the Lyman-a forest 
and the DLA into regions of H I column density which are 
currently poorly constrained (see also Peroux et al. 2003, 
2005). By doing so, we hope to place the SLLS within 
the larger framework of high redshift QSO absorption line 
systems, and to explore their cosmological significance. 
Future papers will examine the ionization state, chem- 
ical abundances, and other physical properties of these 
absorbers. Throughout the paper, we adopt values of the 
cosmological parameters consistent with the latest Wilkin- 
son Microwave Anisotropy Probe (WMAP) results (Ben- 
nett et al. 2003): fl A = 0.7, Q m = 0.3, and H a = 70 km s" 1 
Mpc" 1 . 

2. SPECTROSCOPIC SAMPLE 

The quasar sample in this paper includes spectra from 
two instruments, the Magellan Inamori Kyocera Echclle 
(MIKE; Bernstein et al. 2003) high resolution spectro- 
graph on the Magellan 6.5m telescope at Las Campanas 
Observatory in Chile, and the Echellete Spectrograph and 
Imager (ESI; Sheinis et al. 2002) on the Keck-II 10m tele- 
scope in Hawaii. MIKE is a double echelle spectrograph, 
with a dichroic optical element splitting the beam into 
blue and red arms, each with their own CCD. MIKE pro- 
vides full wavelength coverage, without spectral gaps from 
3350-9500 A in the default configuration. When a 1.0" 
slit is used, MIKE has R = 28,000 and R = 22,000 for 
the blue and red sides, respectively. ESI is a spectrograph 
and imager, which provides continuous wavelength cover- 
age from 3900-10900 A in echellette mode. When a 0.5" 
slit is used, ESI provides R « 9, 000. 

In table 1, we list the 57 QSOs in the current MIKE sam- 
ple. For each QSO in the sample, the data was reduced 
using the MIKE reduction pipeline 8 (Buries, Bernstein & 
Prochaska 2006). The pipeline flat-fields, optimally ex- 
tracts, flux calibrates, and combines exposures to produce 
a single spectrum for the red and blue CCDs of MIKE. In 
Table 2, we list the 56 QSOs in the ESI sample. The bulk 
of these spectra come from the survey for DLA absorption 
presented by Prochaska et al. (2003) but supplemented 
by a new sample of SDSS spectra. All of the ESI data 
were reduced with the ESIRedux pipeline 9 (Prochaska et 
al. 2003). The difference in native resolution of the two 
data-sets will have implications for the Ahi completeness 
limit of each survey. This is discussed in greater detail in 
the following sections. 

2.1. UVES Sample 

In addition to the MIKE and ESI data, we include the 
results of the surveys by Peroux et al. (2003, 2005), which 
we will refer to as the "UVES sample" . All of these data 
were drawn from a heterogeneous sample of high resolu- 
tion observations using the UVES spectrometer (Dckker et 
al. 2000) on the VLT-2 telescope. In five cases, there is an 

8 http://www.lco.cl/lco/magellan/instruments/MIKE/index.html 

9 http://www2.keck.hawaii.edu/inst/esi/ESIRedux/index.html 
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Table 1 
MIKE QUASAR SAMPLE 



Name 


RA (J2000) 


DEC 


(J2000) 




z t 




z 




Q0001-2340 


00 


.03: 


: 15 


.00 


-23 


.2 ! 


: 16 


.5 


2 


262 


1 


.780 


2 


.228 


Q0101-304 


01: 


:03 


:55. 


.30 


-30 


09 


:46 


.0 


3 


.137 


1 


.941 


3 


.095 


SDSS0106+0048 


01 


:06 


:19 


.24 


+00 


:48 


:23 


.3 


4 


.433 


2 


.997 


4 


.378 


SDSS0124+0044 


01 


:24 


:14 


.80 


+00 


:45 


:36 


.2 


3 


807 


2 


.292 


3 


.758 


SDSS0147-1014 


01 


:45: 


:16 


.59 


-09 


45: 


:17 


.3 


2 


.138 


1 


.797 


2 


.106 


SDSS0209-0005 


02 


:09 


:50 


.70 


-00 


:05: 


:06 


4 


2 


.856 


1 


.879 


2 


.816 


SDSS0244-0816 


02 


: 44 


: 17. 


78 


-08 


:16 


:06 


1 


4 


047 


2 


.829 


3 


.996 


HE0340-2612 


03 


:42 


:27. 


.80 


-26 


:02 


: 43 


.0 


3 


082 


2 


.016 


3 


.040 


SDSS0912+0547 


09 


12 


10 


.35 


+05 


:47: 


12 





3 


248 


2 


.146 


3 


.205 


SDSS0942+0422 


09 


:42 


:02 


.04 


+04 


:22 


: 44 


.6 


3 


273 


1 


.896 


3 


.229 


HE0940-1050 


09 


:42 


53 


.10 


-11 


:04 


25 


.0 


3 


.067 


1 


.944 


3 


.025 


SDSS0949+0355 


09 


:49: 


32 


.27 


+03 


35 


31 


7 


4 


.097 


2 


.636 


4 


.045 


SDSS1025+0452 


10 


:25: 


:09 


.64 


+04 


:52: 


:46 


7 


3 


.236 


2 


.108 


3 


.193 


SDSS1032+0541 


10 


:32 


:49. 


.88 


+05 


:41: 


:18 


.3 


2 


.843 


1 


.826 


2 


.804 


CTS0291 


10 


:33 


:59 


.90 


-25 


:14: 


:26 


.7 


2 


552 


1 


.747 


2 


.515 


SDSS1034+0358 


10 


:34 


:56 


31 


+03 


:58 


59 


.3 


3 


.367 


2 


.208 


3 


.322 


Ql 100-264 


11 


:03 


25 


.60 


-26 


:45 


:06 


1 


2 


140 


1 


.747 


2 


.108 


HS1104+0452 


11 


:07 


:08 


.40 


+04 


:36 


:18 


.0 


2 


.660 


1 


.747 


2 


.622 


SDSS1110+0244 


11 


:10 


:08 


.61 


+02 


:44 


:58 


1 


4 


149 


3 


.367 


4. 


.097 


SDSS1155+0530 


11 


55 


:38 


,60 


+05 


:30 


;50 


.5 


3 


.464 


2 


.282 


3 


.418 


SDSS1201+0116 


12 


:01 


: 14. 


.36 


+01 


:16 


:11 


.5 


3 


.215 


2 


.002 


3 


.172 


LB1213+0922 


12 


:15 


:39 


.60 


+09 


:06: 


:08 


.0 


2 


.713 


1 


.796 


2 


.675 


Q1224-0812 


12 


:26 


:37 


.50 


-08 


:29 


:29 


.0 


2 


142 


1 


.747 


2 


.110 


SDSS1249-0159 


12 


:49 


:57 


.24 


-01 


59 


:28 


.8 


3 


662 


2 


.406 


3 


.614 


SDSS1307+0422 


13 


:07 


:56 


.73 


+04 


:22 


:15 


.5 


3 


.026 


1 


.821 


2 


.985 


LBQS1334-0033 


13 


:36 


:46 


.80 


-00 


:48 


:54 


2 


2 


809 


1 


.829 


2 


.770 


SDSS1336-0048 


13: 


;36 


:47 


.14 


-00 


:48: 


:57 


2 


2 


.806 


1 


.796 


2 


.767 


SDSS1339+0548 


13: 


:39 


11 


.95 


+05 


:48: 


:22 


1 


2 


.969 


1 


.969 


2 


.928 


HE1347-2457 


13 


:50 


:38 


.90 


-25 


:12 


:17 


.0 


2 


599 


1 


.749 


2 


.562 


Q1358+1154 


14 


00 


39 


10 


+11 


20 


:22 


.3 


2 


.578 


1 


.747 


2. 


.541 


SDSS1402+0146 


14 


:02 


:48. 


.07 


+01 


:46: 


34 


1 


4 


.187 


2 


.618 


4 


.134 


SDSS1429-0145 


14 


:29 


03 


.03 


-01 


:45: 


19 


.3 


3 


.416 


2 


.331 


3 


.371 


Q1456-1938 


14 


:56 


:49 


.83 


-19 


:38 


:52 


.0 


3 


.163 


1 


.879 


3 


.120 


SDSS1503+0419 


15 


:03 


:28 


88 


+04 


19 


19 


.0 


3 


.666 


3 


.126 


3 


.618 


SDSS1521-0048 


15 


21 


19 


.68 


-00 


:48: 


:18 


.6 


2 


.935 


2 


.178 


2 


.895 


SDSS1558-0031 


15 




10 


.15 


-00 


31 


:20 


.0 


2 


.831 


1 


.784 


2 


.792 


Q1559+0853 


16 


:02 


:22 


.60 


+08 


:45 


:36 


.3 


2. 


.267 


1 


.747 


2 


.218 


SDSS1621-0042 


16 


:21 


:16 


.92 


-00 


:42 


:50 


.8 


3 


.704 


2 


.142 


3 


,656 


Q1720+2501 


17 


:22 


:52 


.90 


+24 


:58 


:34 


.7 


2 


,250 


1 


.961 


2 


,217 


PKS2000-330 


20 


:03 


:24 


.10 


-32 


:51 


:44 


.0 


3 


.776 


2 


.422 


3 


,727 


Q2044-1650 


20 


:47 


:19 


.70 


-16 


:39 


:05 


.8 


1 


,939 


1 


.747 


1 


,909 


Q2053-3546 


20 


:53 


:44 


.60 


-35 


:46 


:52 


.4 


3 


,484 


2 


.154 


3 


.438 


SDSS2100-0641 


21 


:00 


:25 


.03 


-06 


:41 


:46 


.0 


3 


,118 


2 


.093 


3 


.076 


SDSS2123-0050 


21 


:23 


:29 


.46 


-00 


:50 


:52 


.9 


2 


.278 


1 


.837 


2 


.244 


Q2126-158 


21 


:29 


:12 


.20 


-15 


:38 


:41 


.0 


3 


.278 


1 


.943 


3 


.234 


Q2147-0825 


21 


:49 


:48 


.20 


-08 


:11 


:16 


.2 


2 


127 


1 


.879 


2 


.095 


SDSS2159-0021 


21 


:59 


:54 


.45 


-00 


:21 


:50 


.1 


1 


.963 


1 


.797 


1 


.932 


HE2156-4020 


21 


:59 


:54 


.70 


-40 


:05 


:50 


.0 


2 


.530 


1 


.747 


2 


.494 


HE2215-6206 


22 


:18 


:51 


.00 


-61 


:50 


:43 


.0 


3 


.317 


1 


.720 


3 


.273 


Q2249-5037 


22 


:52 


:44 


.00 


-50 


:21: 


:37 


.0 


2 


.870 


1 


.788 


2 


.830 


SDSS2303-0939 


23 


:03 


:01 


.45 


-09 


:39: 


:30 


.7 


3 


.453 


2 


.241 


3 


.407 


HE2314-3405 


23 


:16 


:43 


.20 


-33 


:49 


:12 


.0 


2 


.944 


1 


.684 


2 


.904 


SDSS2346-0016 


23 


:46 


:25 


.67 


-00 


:16 


:00 


.4 


3 


.467 


2 


.166 


3 


.421 


HE2348-1444 


23 


:51 


:29 


.80 


-14 


:27 


:57 


.0 


2 


.933 


1 


.837 


2 


.893 


HE2355-5457 


23 


:58 


:33 


.40 


-54 


:40: 


:42 


.0 


2 


.931 


1 


.854 


2 


.891 



2.187 



2.523 



1.842,2.251 



2.059 



2.279 



"These redshifts correspond to LLS or DLA which were known to exist along the QSO sightlinc prior to the higher resolution observations and 
also which inspired the observations. 
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Table 2 
ESI QUASAR SAMPLE 



NeLIIlG 


RA (J2000) 


DEC 


(J2000) 




z start 


z 


end 


■y' 


i 

m.a.fik. 


Q 1 IQQnm Q i 1 qc:q 
oiJooUUlo-t-loOo 


00 


:13 


:28 


.21 


1 1 q . 

+ 16. 


■58 


:27. 


.0 


3. 


565 


2. 


747 


Q 

o. 


Olo 


o 
a 


-Zo4 


DYflflQ A I 1 R 


00 


:34 


;54 


.80 


i 1 P.. 
+ 16 


39 


:20 


.0 


4. 


290 


3 


.031 


A 

•4. 


onv 
2U i 


A 


ofin 
2oU 


blJbbUUoo+Ullo 


00 


:58 


:14 


.31 


+01: 


15 


:30 


.3 


2. 


535 


2. 


.373 


2 


.499 






biJbbU42 < -UU 


01 


:27 


:00 


.70 


nn. 
— UU. 


,15 


:59 


.0 


4. 


066 


2 


.907 


A 
'1 


m a 

U14 


Q 

o. 


707 
. ( 2 ( 


PQCni QA 1 Q9H7 


01 


:34 


:21 


.60 


1 QQ. 

+oo 


07 


:56 


.0 


4. 


525 


3. 


.154 


A 

•4. 


.4oy 


Q 
O 


7K1 
. / Dl 


blJbbU4oy-Uo24 


01 


:39 


:01 


.40 


no 
— Uo 


24 


:43 


.0 


3. 


008 


2 


.373 


o 

Z. 


yo * 


o 
2 


fi77 
.DM 


bUbbU44Z-t-UU2o 


01 


:42 


:14 


.74 


i nn 
+UU. 


23 


:24 


.0 


3. 


363 


2. 


356 


Q 

o 




Q 

a 


QA 7 
■ o4 / 


blJbbUZ20-t-UU04 


02 


:25 


:54 


.85 


i nn 
+UU. 


54 


:51 


.0 


2. 


963 


2 


.331 


o 
Z. 


noo 


o 
2 


71 /I 
. / 14 


bJJbbUolD+UU4U 


03 


:16 


;09 


.84 


+00: 


40 


:43. 


.2 


2. 


907 


2. 


331 


2. 


.866 






r>l\,JU4Zo-Z2U2 


04 


:26 


:10 


.30 
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"These redshifts correspond to LLS or DLA which were known to exist along the QSO sightlinc prior to the higher resolution observations and 
also which inspired the observations. 
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overlap in quasars observed between our sample and the 
UVES sample. In these cases, we remove those quasars 
which contribute the smaller redshift path to the total 
sample. 

2.2. Redshift Path 

For each quasar in our sample, we define a redshift in- 
terval Az = z en d — z s tart to construct the redshift path 
in a manner similar to that presented by PHW2005. The 
starting redshift z start is given by the lowest redshift at 
which we could identify strong Lya features. For the ma- 
jority of sightlines, this is set by the starting wavelength 
of the spectrum Ao: z sta rt = Ao/1215.67A— 1. Higher 
values for z star t were adopted for a sightline when there 
was either intervening LLS absorption which completely 
removes the QSO flux, or when the SNR of the QSO be- 
came so low as to significantly impact the likelihood of 
detecting a strong Lya feature. The ending redshift z en d 
is given by z enc t — 0.99z em — 0.01, where z em corresponds 
to the QSO redshift. Thus, z en d is defined to be located 
~ 3000 km s _1 blueward of the QSO Lya emission line. 
This offset was chosen to remove those LLS which could 
be associated with the QSO. Tables 1 and 2 list z s t a rt and 
z enc i for each QSO in our sample. Also given in Tables 1 
and 2 are the redshifts, zdla, of any DLA known to be 
present in the data prior to the observations. It is our 
expectation that SLLS are strongly clustered with DLA 
(e.g. Prochaska & Wolfe 1999) and we wish to avoid bi- 
asing the sample (because these systems frequently moti- 
vated the observations). We mask out regions 1500 km s _1 
on either side of the DLA redshifts from the redshift path 
(corresponding to «15 comoving Mpc h~ l at z — 3) to pre- 
vent SLLS clustering with DLA from biasing the sample. 
In this fashion we construct a sensitivity function gi(z) 
which has unit value at redshifts where SLLS could be de- 
tected and zero otherwise. The combined redshift path of 
our sample is Az = 124.4. 

The H I frequency distribution function /hi (A, A) for 
the SLLS describes the number of SLLS in a range of col- 
umn densities (N,N + dN), and a range of absorption 
distances (A, X + dX), 

f m (N,X)dNdX, (1) 

where the absorption distance (Bahcall & Peebles 1969) is 
given by 



AX = 




By considering dX instead of dz, /hi (AT, X) is defined over 
a constant, comoving pathlength but we have introduced 
our assumed cosmology into the analysis (e.g. Lanzetta 
1993). The cosmological term in fm(N,X) is still rel- 
atively obscure, but can be made less so when we con- 
sider that for z > 2, w + z) /f2 m for a standard 
ACDM model. By way of example, consider a sightline 
with Az — 1 pathlength at the average redshift for our 
survey z — 3.1 corresponding to AX 3.7. 

To characterize the survey size, we construct the to- 
tal redshift sensitivity function, g(z) — E^i(z). Figure 1 
shows the g(z) curves for the MIKE, ESI and UVES sam- 
ples, and the combined g(z) for all three samples. We find 



that the total sample presented here increases g{z) by fac- 
tors of 2 to 8 over the redshift range 2 < z < 4 when com- 
pared to the previous survey of Peroux et al. (2005) when 
we consider the column density range 19.3 cm _2 <log Njji 
< 20.2 cm~ 2 . In Figure 1, there is a sharp feature in g(z) 
at a redshift of z ~ 2.7. This feature is due to a gap in 
the ESI data at A ~ 4500A which appears in every ESI 
spectrum in the sample. At the redshifts surveyed, the 
resolution of the data is sufficiently high that sky lines do 
not effect our ability to detect LLS in the data. 

3. A H i analysis 

With the survey size defined, we now turn to the SLLS 
in the survey and discuss the measurement of the Ahi 
for the SLLS, and an estimate of the completeness of the 
sample. 

3.1. Identifying SLLS and Measuring Nja 

To identify the SLLS in our survey, every spectrum from 
the ESI and MIKE samples was visually inspected for 
strong absorption in H I Lya. When a candidate SLLS 
absorber is identified, we first assign a local continuum 
level to regions approximately 2000 km s _1 on either side 
of the absorption. We do not first continuum normalize 
the data, opting instead to model the absorption and con- 
tinuum simultaneously. This is done because the damping 
wings of the Lya profile may depress the continuum for 
several tens of Angstroms. It is also necessary to vary the 
continuum level along with a model to estimate the errors 
on the Nm values. When the absorption occurs near the 
QSO emission line features and/or when the redshift of 
the absorption is high, the assignment of the continuum 
level is often the dominant source of error. In this work, 
we assign a minimum value of the error on the log Nhi to 
be 0.05 dex from continuum level errors alone. 

Next, we compared a Voigt profile with log Njji = 19.0 
cm~ 2 at the center of the absorption to determine if it 
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Fig. 1. — Redshift sensitivity function g(z) for the samples. The 
sharp feature at z ~ 2.7 is caused by a gap in the ESI data. The 
upper (black) curve represents the full sample g(z)subject to the 
added constraint that it apply to the data where a SLLS with Nm 
> 19.3 cm~ 2 could be found, with the other curves (blue, red, green) 
representing the g(z) for the MIKE, ESI, and UVES data samples 
which comprise the full sample. 
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Fig. 2. — H I Lyman-a transitions for the SLLS in the MIKE sample. The velocity zeropoint is determined from low ionization metal lines 
whenever possible. Overlaid (in red) are the single component Voigt profile models for the absorption along with the ±1<7 error estimates (in 
green) on the Nm. The horizontal dashed line represents the local continuum level. The vertical dashed line highlights the velocity zeropoint. 
For each absorber , the QSO name, absorption redshift, and Nm estimate are listed at the top of each panel. In the case of blended SLLS 
where the individual Lyc* lines are separated by more than 300 km s" 1 , we show the individual fits (e.g. the SLLS in PKS2000-330). 
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Fig. 2 - continued 
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was consistent with at least this amount of atomic hy- 
drogen gas. If metal-line transitions were identified, we 
use the redshift of the metals to more accurately assign 
a redshift for the gas. We do not, however, require that 
metal line absorption be present, because this would in- 
troduce an undesirable bias. For most of the absorption 
systems in our sample, the data cover a wide variety of 
metal line transitions for each absorber. Common transi- 
tions include those of C II, C IV, O I, Si II, Si IV, Al II, 
and Fc II. In the cases where both low and high ionization 
metal line transitions are present, we choose to use the red- 
shift given by the low-ion absorption, because these ions 
more likely trace the atomic hydrogen gas. Finally, for 
those absorbers meeting our minimum Ahi condition, we 
fit the absorption by stepping through values of Ahi un- 
til the absorption is well modeled. The Lya line profile 
and continuum level were modeled using custom software 
(the X-fitdla tool within the XIDL package 10 ). In the cases 
when two or more SLLS occur within Sv < 300 km s _1 of 
each other, the total Ahi is reported instead of individual 
Ahi measurements. This is because the individual Ahi 
measurements are highly degenerate and because the gas 
may physically bound to a single virialized halo. 

For each absorber, a number of effects contribute to the 
error in the Ahi value. As mentioned above, we believe 
that errors in the continuum level contribute at least 0.05 
dex to the error. The continuum level error tends to in- 
crease with redshift because the amount of absorption from 
Lyman-a forest increases and one is less certain to identify 
regions free of absorption. The increase in Lyman-a for- 
est lines also affects the Am measurement in the damping 
wings and line core regions of the Lya line because of en- 
hanced line-blending. Redshift uncertainties were only a 
minor source of error for each absorber, since nearly every 
absorption system showed at least one metal-line transi- 
tion. Finally, the Poisson noise of the data adds additional 
uncertainty to the fitted value. 

In Figures 2 and 3, we show the profile fits and the ±lcr 
fits for each LLS in our sample. The values for the z, Ahi, 
and ct(Ahi) are given in Tables 3 and 4. The LLS listed 
in these tables do not represent every LLS present in our 
data. In a number of cases, we removed an SLLS from the 
sample because the QSO was specifically targeted to study 
this SLLS. For those SLLS removed, we also remove 1500 
km s _1 of pathlength on cither side of the SLLS in the 
same manner as for the DLA. In total, we have a homoge- 
neous sample of 47 SLLS, of which 17 come from the ESI 
sample and 30 come from the MIKE sample. Including 
the UVES sample, there are now a total of 78 SLLS. In 
Figure 4, we show a histogram of the Ahi values for the 
individual and complete samples. 

3.2. Completeness 

To gauge the completeness of our Ahi analysis we per- 
formed a number of tests. Our primary completeness con- 
cern is with the detection efficiency at the low H I column 
density limit of the ESI sample, where the data has lower 
spectral resolution. At lower resolution, the effects of line 
blending and continuum level placement tend to wash out 
the damping wings of log Nhi — 19.0 Lya lines. As the col- 
umn density increases, the equivalent width of the Lya line 

10 http://www.ucolick.org/~xavier/IDL/ 



Table 3 

LLS FROM THE MIKE SAMPLE 



Name 




abs 


log. 


Vhi 


a(N m ) 


Q0101-304 


2. 


908 


19. 


10 


0. 


.15 


SDSS0106+0048 


3 


286 


19. 


05 


o 


.25 


SDSS0106+0048 


3. 


321 


19. 


10 





20 


SDSS0106+0048 


4. 


172 


19. 


05 





.20 


SDSS0124+0044 


3 


078 


20. 


20 





.20 


SDSS0912+0547 


2. 


522 


19. 


35 


o 


.20 


SDSS0949+0355 


3. 


311 


19. 


75 


0. 


.15 


SDSS1034+0358 


2. 


849 


19. 


.50 


o 


.25 


SDSS1034+0358 


3 


.059 


19. 


15 


o 


.15 


SDSS1034+0358 


3. 


.003 


19. 


10 


o 


.15 


Q1100-264 


1 


.839 


19. 


40 


o 


.15 


HS1 104+0452 


2. 


601 


19. 


.90 





20 


LB1213+0922 


2. 


.523 


20. 


20 


o 


.20 


SDSS1307+0422 


2. 


250 


20. 


00 


o 


.15 


SDSS1402+0146 


3. 


456 


19. 


.20 


o 


.30 


Q1456-1938 


2. 


351 


19. 


55 


o 


15 


Q1456-1938 


2. 


169 


19. 


75 





20 


SDSS1558-0031 


2. 


630 


19. 


50 





.20 


SDSS1621-0042 


3. 


104 


19. 


70 





.20 


PKS2000-330 


3 


188 


19. 


.80 





.25 


PKS2000-330 


3 


172 


19. 


80 





.15 


PKS2000-330 


3 


192 


19. 


.20 





.25 


Q2053-3546 


2. 


350 


19. 


60 





.25 


Q2053-3546 


2. 


989 


20. 


10 


0. 


.15 


Q2053-3546 


3 


094 


19. 


00 


0. 


.15 


Q2053-3546 


2. 


333 


19. 


30 


0. 


.25 


Q2126-158 


2. 


638 


19. 


.25 





.15 


Q2126-158 


2. 


769 


19. 


20 





.15 


HE2314-3405 


2. 


386 


19. 


00 


0. 


20 



Table 4 

LLS FROM THE ESI SAMPLE 



Name 


z 


abs 


log 7Vjji 


a(N m ) 


PX0034+16 


3 


754 


20 


.05 


0. 


.20 


SDSS0127-00 


2. 


.944 


19 


.80 





15 


PSS0808+52 


3. 


.524 


19 


.40 


0. 


.20 


SDSS0810+4603 


3. 


472 


19 


.90 





.20 


SDSS0826+5152 


2. 


862 


20 


.00 





.15 


SDSS1004+0018 


2. 


746 


19 


.80 





.20 


PSS1248+31 


4. 


075 


19 


.95 





15 


PSS1253-02 


3. 


603 


19 


.50 





.15 


SDSS1257-0111 


2. 


918 


19 


.95 





.15 


Q1337+11 


2. 


508 


20 


.15 





.15 


PSS1432+39 


3 


994 


19 


.60 





.25 


FJ2129+00 


2. 


735 


20 


.10 





20 


PSS2241+1352 


3. 


654 


20 


.20 





.20 


SDSS2315+1456 


3. 


135 


19 


.95 





15 


SDSS2315+1456 


2. 


.943 


19 


.35 





.20 


PSS2323+2758 


3. 


267 


19 


.40 





.20 


PSS2323+2758 


3. 


565 


19 


.30 


0. 


.20 



also increases and one generally derives a more reliable Ahi 
measurement. Because the MIKE and UVES samples are 
of higher resolution, the effect at lower column densities 
is minimized. After experimenting with mock ESI spectra 
(see below), we chose to limit the ESI analysis to those 
absorbers with 19.3 cm _2 < log Af//< 20.3 cm~ 2 . For the 
MIKE sample, we consider the same range of column den- 
sities as in the UVES sample, 19.0 cm" 2 < log Nhi < 20.3 
cm~ 2 . Our tests indicate that we have detected all SLLS 
in these ranges in Ahi at a 99% level of completeness. 

To verify that the lower bound of log Nhi= 19.3 cm~ 2 
is appropriate for the ESI sample, we randomly placed 
mock LLS onto 50 of the ESI spectra. The column den- 
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Fig. 3. — Same as for Figure 2, but for the SLLS in the ESI sample. 
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Fig. 3 - continued 



sities of these LLS were chosen at random to lie in the 
range 18.0 <logA#/< 20.5 cm~ 2 . In all cases, SLLS with 
log Nhj> 19.3 cm -2 were identified. Furthermore, the col- 
umn density fit to these SLLS were within the 2— a error 
estimates in all cases but one out of fifty. We have also 
used simulated LLS to address the issue that blending of 
multiple LLS with log Nhi < 19.3 cm -2 could appear as a 
single absorber with log Nhi > 19.3 cm~ 2 in the ESI sam- 
ple. As before, we placed mock LLS upon ESI spectra, but 
in pairs with a random velocity separation 6v subject to 
the constraint Sv < 300 km s _1 . This constraint matches 
our velocity separation constraint for the observed SLLS 
as to whether the Ahi is reported as a sum, or as separate 
systems. The blended, mock LLS absorption was then fit 
with a single absorption profile, and the Ahi compared to 
the sum of the individual component Ahi values. We find 
that a pair of LLS with log Nhi< 18.8 cm -2 or less can not 
blend in a manner so as to be fit as an SLLS with log Nhi 

> 19.3 cm -2 . However, a blend of a SLLS with log Ay/ 

> 19.3 cm _2 with a LLS with log Nhi ~ 19.1 cm~ 2 gener- 
ally leads to a systematic overestimate of Ahi- Similarly, 
a blend of LLS with log Nhi= 19.1 cm~ 2 and log Nhi = 
18.8cm~ 2 may mimic an SLLS with log Nhi = 19.3 cm~ 2 . 
In essence, this effect leads to a Malmquist bias for the 
sample at low Ahi value. It is not generally possible to 
directly explore these final two issues without good knowl- 
edge of the incidence frequency, clustering properties, and 
column density distribution of the LLS. Nevertheless, we 
have some constraints from both the MIKE and UVES 
samples. Because these samples are of higher resolution, 
we can often distinguish blended absorption of LLS from a 
single LLS. Moreover, if two LLS are blended, their metal 
lines should give some clues that there is blending, pro- 
vided both LLS have sufficient metallicity. To the extent 
that we can currently distinguish between blended and un- 
blended LLS, we do not believe that the effect of blended 
LLS strongly effect our measurements of the Nhi but we 
caution that a larger sample than the one presented here 
is required to accurately address this issue. 

4. RESULTS AND DISCUSSION 

In this section we discuss a number of results related to 
the H I frequency distribution of the SLLS. For the major- 
ity of the section, we consider two Ahi groups driven by 
the completeness limits of the spectra: (i) 10 190 cm -2 < 
Am < 10 20 3 cm -2 which does not include the ESI sample; 
and (ii) 10 193 cm~ 2 < Ahi < 10 20 3 cm~ 2 which includes 
all of the samples. The groups cover an integrated absorp- 
tion pathlength of AA 19 - = 329.1 and AA 193 = 467.7 
respectively. We maintain this division, as opposed to 
combining the low Ahi results from the echelle data with 



the Nhi > 10 cm analysis, for the following reasons: 
First, the redshift distribution of the ESI survey (Figure 1) 
is significantly higher than that of the UVES data and es- 
pecially the MIKE sample. If there is redshift evolution 
in /hi (AT, X), then it would be erroneous to mix the two 
groups. Indeed, PHW05 find significant evolution int he 
normalization of /hi(A, X) for the damped Lya systems 
and one may expect similar evolution in the SLLS popu- 
lation. Second, the division allows us to investigate sys- 
tematic errors between the various surveys, including the 
UVES analysis. Finally, we will find it instructive valuable 
to consider these two Ahi groups separately when focusing 
on the behavior of /hi (A, X) near Ahi = 10 19 cm" 2 . 

4.1. Power-Law Fits to the /hi(A, X) Distribution of the 

SLLS 

In Figures 5 and 6, we show the binned evaluation of 
/hi (A, X) for each of the data-sets, and their combined re- 
sults. The results are summarized in Table 5. It is evident 
from the figures that the /hi(A, A) distributions can be 
reasonably modeled by a power-law over the SLLS regime: 



/hi (A, A) = k 1 N° 



(3) 



We have performed a maximum likelihood analysis on the 
unbinncd Ahi distribution to constrain a and set the nor- 



-1.40 



+0.15 



and 



malization of /hi (A, A). We find aigo — -i.™_ 16 
ai93 = — 1.191q' 2 ° for the combined results of the two 
groups. Note that the reported errors do not include co- 
variance terms. The results for the individual data samples 
are all consistent (within 2a) with the combined results 
and one another (Table 5). 




19.0 



19.5 



20.0 



l°g N HI 

Fig. 4. — Histograms showing the number of systems detected as 
a function of log Njji for the various samples. The upper two panels 
represent the MIKE and ESI data, which are presented for the first 
time in this paper. The third panel shows the data from Peroux et 
at. (2005). The bottom panel shows the full Nm histogram when 
we combine the three data sets. 
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Table 5 
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^■ 48 °-0.062 
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-21.48+°;i| 


0.182 


+0.027 
-0.023 
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4^ 
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a = -1.40 
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logN ra 
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Fig. 5. — Values for the H I frequency distribution for the SLLS 
with 19.0 < logN HI 20.3cm -2 . Shown separately are the MIKE 
data presented first here, the UVES sample of Peroux et al. (2003, 
2005), and the two samples combined. The four bins span equal size 
Alog Nh i= 0.3 except for the last bin which has Alog Nh i— 0.4. 
Also shown is the result of a maximum likelihood analysis on the 
unbinned data for a power-law model of fm(X, X). The power-law 
index a is given for each sample. 
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Fig. 6. — Same as for Figure 5, but with the inclusion of the ESI 
data, and with the range 19.0 < log Nfji 20.3 cm -2 , over 3 bins 
equal in size to the 3 highest H I column density bins in Figure 5. 
Also shown is the result of a maximum likelihood analysis on the 
unbinned data for a power-law model of fm{X, X). The power-law 
index a is given for each sample. 
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One notes that the power-law is shallower for the AT^- 3 
sample than the sample. Although the effect is not 

statistically significant, we do find the trend is apparent in 
the independent MIKE and UVES samples. These results 
suggest that the fm(N,X) distribution is steepening at 
Ahi < 10 193 cm~ 2 . We will return to this point in § 4.4. 

4.2. Redshift Evolution 

Over the redshift interval z = 2.2 to 3.5, PHW2005 
found that the shape of /hi (A, X) for DLA is roughly con- 
stant but that the normalization (parameterized by the 
zeroth moment, ^dla) increases by a factor of approxi- 
mately two. Similarly, studies of the Lya forest indicate 



(4) 



4llsW= / fm(N,X)dN 

NsLLS 

with N SLLS = 10 190 or 10 19 3 cm- 2 . We calculate 
WPQ(i90) = 0.134l[i; 9 !and^LLSpO(i93) = 0.110±g;g? 



for z < 3and£ SLLS (X) 



(190) 



0.254; 



-0.050 
-0.042 



and e ShhS (X) 



(193) 



0.1261^;^ for z > 3. As is evident from Figure 7 (i.e. fo- 
cus on the values for the lowest TVjji bins), there is no 
significant redshift evolution in ^slls(A') for the Ahi = 
cm~ 2 sample. There is, however, an indication of 



19.3 
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increasing £sllsPQ with increasing redshift for the Ahi 
(z) Lya oc (1 + z) 2 - 6 (Cristiani et al. 1995) at z w 3 imply- 10 19 cm -2 sample. This result is driven by the lowest Am 



ing £(X) Lya oc (1 + z) 21 . One may expect, therefore, that 
^SLls(A') w iH also increase with redshift over this redshift 
interval. 

In Figure 7, we show the results for fm(N, X) in the two 
A/hi groups, but broken into the redshift regions above 
(thin line) and below (thick line) z = 3. A visual in- 
spection of Figure 7 reveals little explicit difference in the 
shape of fnj(N,X) in the two redshift regimes. More for- 
mally, a KS test for the two groups returns probabilities of 
-P190 = 0.189 and P193 = 0.809 that the two redshift sam- 
ples are drawn from the same parent population. There- 
fore, we contend there is no large redshift evolution in the 
shape of fm(N, X) for the SLLS over the redshifts consid- 
ered here. 

Granted that the shape of /hi (N, X) for the SLLS is not 
steeply evolving, we can examine redshift variations in the 
normalization by examining the zeroth moment 
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Fig. 7. — Values for the H I frequency distribution for the 
SLLS separated into two redshift samples, with the thin symbols 
corresponding to z < 3 and the thick symbols corresponding to 
z > 3. The upper panel shows the ESI+MIKE+UVES sample and 
the lower panel the MIKE+UVES sample. The lack of significant 
change in shape of /hi^i^O for the two redshift bins indicates a 
lack of strong evolution with redshift in fm(N, X). 



bin which we caution is the most sensitive to the effects of 
line-blending and that such effects are heightened at higher 
redshift. In summary, there is only tentative evidence for 
redshift evolution in ^sllsO^O- As noted above, this con- 
tradicts the apparent evolution in £(X) for the DLA pop- 
ulation and the Lya forest. If confirmed by future studies, 
perhaps modest evolution in the SLLS population suggests 
the competing effects of the decrease in the mean density 
of the universe and a decrease in the intensity of the ex- 
tragalactic UV background radiation field. 

4.3. Is f m (N,X) for the SLLS flatter than the DLA? 

There is a significant mismatch in power-law exponents 
for the Lya forest (ctLya — —1-5) and the DLA (oldlA ~ 
—2) at z « 3. This difference predicts that the shape of 
fm{N,X) change appreciably between these two column 
densities, i.e., the fm(N, X) distribution of SLLS is flatter 
than that of the DLA. 

The simplest test of this prediction is to measure the 
power-law index of /hi (AT, X) in the SLLS regime and com- 
pare against its slope near A^hi = 10 20 ' 3 cm~ 2 . Both the 
Gamma function and double power-law fits to fm(N,X) 
for the SDSS DLA sample indicate udla = —1-9 ± 0.1 at 
AT H i < 10 21 cm" 2 (PHW2005). We compare these values 
with asLLS = -1-19 ± 0.21 for log N l ™ = 19.3 cm" 2 and 
asLLS = -1.40 ± 0.15 for log ATgp = 19.0 cm- 2 . The 
statistical difference in the power-law slope is greater than 
2a significance but not a full 3c result, in terms of the 
exponents alone. 

Figure 8 shows the /hi (AT, X) distributions for the SLLS 
and SDSS DLA to a column density of log Nhi— 21.5cm- 2 . 
As before, we show the results for the different SLLS AThi 
groups separately. Overplotted on the data are the SLLS 
single power-law fit and the low column density end of 
the double power-law for the SDSS-DLA (PHW2005). It 
is apparent that a simple extrapolation of the SDSS-DLA 
fit significantly overpredicts the frequency of the SLLS, 
especially at \ogNHi< 19.7 cm -2 . We contend, there- 
fore, that the /hi (AT, X) flattens in slope around the Ahi 
corresponding to the canonical definition for the DLA, 
log N HI = 20.3cm- 2 . Zheng & Miralda-Escude (2002) pre- 
dicted a flattening in /hi (A", X) based on photoionization 
and self-shielding models of isothermal gas profiles in dark 
matter halos. Although their analysis examined individ- 
ual halos at a given mass, Zheng (2006, priv. comm.) has 
convolved their results with a Sheth-Tormen halo mass 
function (Sheth & Tormen 1999) and predict f m (N,X) 
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should flatten at Am w 10 cm , consistent with our 
results. 

There are a few cautionary remarks to make regard- 
ing Figure 8. First, the shape of the selection functions 
g(z) for the SLLS and DLA samples do not exactly match 
because the SLLS database includes a somewhat higher 
fraction of z > 3 quasars, and thus a fractionally larger 
g(z) at those redshifts. Although the mean differs by only 
Sz — 0.3, we note that the comparison is not perfect as 
we have not considered any evolution in the normaliza- 
tion of fm(N,X) for the SLLS or DLA but have simply 
plotted the full samples. Another systematic effect is that 
the SDSS-DLA sample may suffer from a Malmquist bias. 
Specifically, the statistical and systematic errors (e.g. the 
effects of line blending) in the Am values of the DLA are 
significant and will drive the observed /m(A, X) distribu- 
tion to a steeper slope. It is possible that this effect ex- 
plains the marked drop in /hi (A, X) at Ahi = 10 20 7 cm~ 2 
in Figure 8. We intend to address this issue directly with 
follow-up, higher resolution observations of a large sam- 
ple of SDSS-DLA with N m « 10 2a3 c m - 2 . If there is a 
substantial Malmquist bias in the DLA sample, then the 
decrease in the slope of fm(N,X) would be more gradual 
than that suggested by Figure 8. 

4.4. Is There an Inflection in fm(N,X) within the SLLS 

Range? 

While the mismatch between the DLA and the Lya for- 
est in the power-law description of their fm{N, X) dis- 
tributions suggests that /hi (A, X) for the LLS will show 
intermediate values (a « —1-7), the observed incidence 
of LLS reveals a different result. As PHW2005 discussed, 
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Fig. 8. — Values for the H I frequency distribution for both the 
SLLS and the DLA sample of PHW2005. Overplotted are the results 
for the single power-law for the SLLS from this work (dot-dashed 
line, power-law index a = —1.19 and a = —1.40 for the upper and 
lower plots, respectively), and the low column density end of the 
double power-law fits to fni(N,X) for the DLA from PHW2005 
(dashed line, power law index alpha = —2.0). Neither fit describes 
the full range well, and the fm(N,X) shows a flattening near the 
canonical DLA definition of logN[ji= 20.3cm~ 2 . 



a simple spline interpolation of the DLA and Lya forest 
/m(A, X) distributions through the LLS regimes predicts 
over an order of magnitude more LLS than observed per 
AX. PHW2005 argued, therefore, that the f m (N,X) dis- 
tribution for the LLS must exhibit an inflection as ev- 
idenced by dlog //g? log A > —1.5. Zheng & Miralda- 
Escude (2002) have also predicted that there should be 
an inflection in /hi (A, X) in the SLLS regime for galaxies 
exposed to an ionizing radiation field. It is worth investi- 
gating with our data-set whether evidence exists for just 
such an inflection. 

The simplest approach is to examine whether d log f/d log N > 
— 1.5 in the LLS regime. Regarding our results on the 
SLLS, we find that dlogf/dlogN > -1.5 for both the 
logAgf = 19.0 cm" 2 and logX^f = 19.3 cm" 2 groups 
(Table 5). The differences, however, have less than 2er 
significance. Using only the current data-set and the dis- 
tribution of DLAs, we do not report the existence of an 
inflection in /m(X, X) within the SLLS regime. 

4.4.1. Constraints from lower Am LLS 

In order to pursue the question further, we introduce two 
new observational constraints on the LLS. The statistical 
significance for the SLLS alone is limited by the combi- 
nation of sample size and observed baseline in log Nhi- 
We cannot arbitrarily increase the sample size, but we are 
able to introduce new constraints which are sensitive to 
lower column density LLS. The number density of opti- 
cally thick LLS has been well constrained by many stud- 
ies. For our purposes here, we apply the constraint at 
redshift z = 3 of £(X) LLS = 0.7 ± 0.1 (Sargent, Steidel 
& Boksenberg 1989; Storrie-Lombardi et al. 1994; Peroux 
et al. 2003). We also use a measure of the incidence of 
optically thin partial Lyman Limit Systems (PLLS) from 
Buries (1997), who found 12 systems with mean redshift 
at z — 3 with 17.2 cm~ 2 < logXjy/ < 17.8 cm" 2 over a 
redshift path of Azplls = 59.07 and absorption path of 
AXp LLS = 16.63. This gives £(X)plls — 0.20 ±0.06. 

4.4.2. The observable distribution function 0(X, X) 

In order to facilitate the analysis and interpretation of 
the distribution over the full range of LLS and DLAs, 
we introduce a function called the observable distribution 
function of H I: 

m t{X) 



G{N,X) 



(5) 



A log A AX A log TV ' 
where m is the number of systems observed over an ab- 
sorption path, AX, and a column density range A log A. 
This is simply the frequency distribution function in log- 
arithmic column density bins and is related to the classic 
distribution by: 

0(N, X) dlogN = /(A, X)dN , 0{N, X) = f(N, X) N ln(10) 

(6) 

The observable distribution has a few nice features. It is 
unitless, and gives the direct number of systems observed 
over a specified bin in logarithmic column density (almost 
all studies of H I absorption show and analyze the data 
in bins of constant width in log Nhi)' It a lso removes 
one factor of column density from the steep slope in the 
frequency distribution, which enables better assessments 
of change in slopes, as well as smaller effects from the 
rapid change in the distribution over bins of large size. 
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4.4.3. C(N,X)for the LLS and the DLA 

In Figure 9, we show the observational constraints on 
0(N, X) from the present sample of SLLS and the DLA 
from PHW2005. In addition, we show an observational 
constraint on the abundance of partial Lyman limits at 
z = 3, Oplls — 0-42 ± 0.13, corresponding to the column 
density interval (17.2 cm~ 2 < log N HI < 17.8 cm~ 2 ; Buries 
1997). 

We fit a simple analytic function to the combined set 
of ALL-19.0 SLLS (Table 5), the PHW2005 DLA sample 
for logA^H7> 20.3 cm -2 , as well as the two additional 
LLS constraints described above. We fit a three parameter 
model expressed as: 

log 0(N, X) = a + b* (log N HI - 19) - 10 log Nhi ~ c (7) 

This parameterized form can be recast into a more familiar 
r-distribution (e.g. Pei & Fall 1995; Peroux et al. 2003): 
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Fig. 9. — a) Values for the observable distribution function for the 
entire range of H I absorbers which are optically thick at the Lyman 
limit. The solid line is a best fit 3 parameter function (assuming 
Gaussian statistics) of the data which include the SLLS presented 
here, the DLA from PHW2005, and a bin which contains the con- 
straint on the number of partial LLS. The dot-dashed line shows 
the best fit after removing the two outlying bins from the PHW 
sample. Also shown (dotted line) is the 0(N, X) for the Lyman-o 
forest, b) Same as a), but extended to include the Lyman-a forest 
results of (Kirkman &: Tytler 1997), which are shown as the trian- 
gles. The constraint from the number of partial LLS is displayed by 
the square. Overlaid are the values for the power-law slope a for 
the TVjji regions: ctLya corresponding to the optically thin Lyman— 
a forest with 12.1 < log Nhi < 14, aoTLya corresponding to the 
optically thick Lyman-a forest with 14 < log Nhi < 16, culls 
corresponding to 16 < log Nhi< 20.3, and oldla corresponding 
to logTVjji = 20.3 — 21.3, to highlight the observed and predicted 
changes in the logarithmic slope of fm(N, X). 



0(N,X) = O a ^y +1 (8) 

where log AT, = c - 0.3622, logO = 6 log A* + (a - 
196), a = 6-1. 

We show the best fit function to the binned data as- 
suming Gaussian statistics as the solid line in Figure 9, 
parameterized by (a, 6, c) = (—0.621, —0.189, 21.51), which 
gives a reduced % 2 for 17 degrees of freedom (dof), xl = 
1.86. We assess the dependence of the fitting on the in- 
cluded data-sets, by sequentially removing constraints. If 
we do not include the constraint on the total number 
of optically thick Lyman limits, the best fit parameters 
are virtually unchanged, and the fitted parameters above 
give a predicted number of optically thick Lyman limits 
of: £(X)lls = 0.68. The two outliers in the above fit 
are the DLA points centered near bins of log Nhi— 20.3 
and 20.6 cm -2 . If we drop the first of these points, we 
find (-0.660,-0.198,21.56) and xl = 1-24 for 15 dof. 
Dropping both DLA bins gives an acceptable fit, with 
(-0.677,-0.227,21.63) and xl = 0.552 for 14 dof. The 
last fit is shown as a dashed-dotted line in Figure 9. 

The fits show that the 0(N, X) slope in the SLLS region 
falls between -0.227 < b < -0.189 (recall a = b - 1). We 
also find that the HI cutoff scale is between 21.51 < c < 
21.63 which matches the results of PHW2005 and Peroux 
et al. (2003). All the fits produce a reasonable number 
of optically thick Lyman limit systems: 0.61 < ^lls(A) < 
0.68. 

4.4.4. Are the SLLS a distinct population? 

Finally, we present another moment of the Am distribu- 
tion, NmO(N, X), in Figure 10. This represents the total 
H I column density per unit logarithmic column density 
per unit absorption path. The SLLS and DLA data points 
and the first analytic fit (to the entire data set) is shown. 
The data are relatively flat from Am = 10 20 3 cm~ 2 to 
10 215 cm' 2 with the functional form peaking at logNm = 
21.1. It is evident that absorbers with Am ~ 10 21 cm~ 2 
dominate the mass density of H I in the universe. Con- 
sider a comparison of the SLLS and DLA. Whereas the 
damped Lya systems with 20.3 < log Am < 21.5 domi- 
nate the mass density with roughly equal contribution per 
A log N, the SLLS lie distinctly below, with the absorbers 
at Am = 10 19 cm~ 2 adding a negligible portion and the 
full SLLS range in Am contributing only w 15% to the 
total mass density. We note that Zwaan et al. (2005) ob- 
tain a very similar shape for AhiO(A, A)at z = 0. Again, 
this behavior is likely related to the fact that the majority 
of SLLS are highly ionized. The results in Figure 10 lends 
further support to the concept that the SLLS absorbers 
are a distinct population from the damped Lya systems. 

4.5. Implications for lower Ahi 

Armed with a description for the full range of H I col- 
umn densities which are optically thick at the Lyman limit, 
we now wish to extend our analysis downward in Am to 
the Lyman-a forest. Although a more thorough analy- 
sis is required to fully assess the statistically acceptable 
distributions over such a large range in log Nhi, we can 
clearly show that slopes in the SLLS region as steep as 
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a = — 1.5 are unacceptable. In Figure 9b, we overplot the 
extrapolation of the best fit observable distribution of low 
column density Lyman-a absorbers (12.1 < log Nhi < 14) 
as a dotted line with a slope of b = —0.5 (i.e., fm(N, X)oc 
TV- 1 - 5 ; Kirkman & Tytler 1997). Although the extrapo- 
lation has uncertainties related to the normalization and 
completeness, it does highlight the overprediction of SLLS 
based on a simple power-law extrapolation from lower col- 
umn density studies. Even if we allowed for freedom in 
normalization, it is clear that in no region except near 
log Nhi— 20.5 cm~ 2 is a = —1.5 a good description of the 
high log Nhi distribution. A subset of the data which were 
used constrain the fit to Lyman-a absorbers is shown in 
the lower panel of Figure 9 as triangles. 11 A gamma func- 
tion with low-end slope of a = —1.2 is a good fit to the 
high density data sets of LLS and DLAs. But the extrap- 
olation of the low end slope to the regime of the Lyman-a 
absorbers underpredicts the observed numbers by almost 
a factor of 100! 

In order to reconcile the data presented here, together 
with the DLAs, LLSs and the Lyman-a absorbers, the full 
distribution function must contain at least three changes 
in logarithmic slope of the frequency distribution (or in- 
flections, d 2 f/dN 2 = 0), as previously argued (in part) 
by Bechtold (1987) and Petitjean et al. (1993). The 
first inflection is seen in the change between the DLA 
distribution and SLLS. This is required if one includes 
constraints from the LLS and the PLLS and demands 
constant slopes down to these column densities. There 
must then be at least two more changes in the logarith- 
mic slopes, one to account for the much higher number 
of Lyman-a absorbers, and the second to finally merge 
back to the Lyman-a slope with a — —1.5. We present 
one such solution with a dotted line bridging the gap be- 
tween log Nhi= 14 cm -2 and log Nhi— 16 cm~ 2 with a 

11 Note, the two points above log Njji= 14.5 cm~ 2 were not included 
in the power-law fit. 




19.0 19.5 20.0 20.5 21.0 21.5 22.0 
log N HI 



Fig. 10.— Values for N m O(N,X), the first moment of the ob- 
servable distribution function over the SLLS and DLA range. The 
data give the relative contributions to the cosmic density of hydro- 
gen atoms per unit log Nhi at z ~ 3, and show that DLA with 
log N[ij= 20.3 — 21.5 overwhelmingly dominate the H I mass den- 
sity, with the SLLS contributing only 15% to the total. 



frequency slope of a = —1.9. The Lyman-a absorbers are 
claimed to have a single power law slope over two decades 
of A r Hi(12.5cm _2 < logNni <14.5cm~ 2 ), and are classi- 
fied as a single population. In comparison, the absorbers 
spanning 17cm~ 2 < logNni <20cm~ 2 , over three decades 
of AThi, are well described by a single power law slope of 
a = —1.2 and also could be classified as a single population 
by applying the same argument. 



5. concluding remarks 

The SLLS results presented in this paper motivate sev- 
eral avenues of future exploration. First, we must extend 
the survey downward in LLS H I column density. An im- 
proved constraint on the number of optically thick LLSs 
would significantly tighten the constraint on the low-end 
slope of the SLLS distribution. Second, a large sample of 
higher H I column density Lyman-a forest data is needed. 
Specifically, a large sample of Lyman-a forest with 14.5 < 
logNni < 16.5 cm~ 2 directly tests the predicted shape 
of the H Idistribution. Finally, one must pursue detailed 
ionization studies of the LLS to fully assess their baryonic 
mass, metallicity, etc. Our survey will provide the data-set 
required for just such studies. 
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